Effect of disorder close to the superfluid transition in a two-dimensional Bose gas 
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We experimentally study the effect of disorder on trapped quasi two-dimensional (2D) ^'^Rb 
clouds in the vicinity of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition. The disorder 
correlation length is of the order of the Bose gas characteristic length scales (thermal de Broglie 
wavelength, healing length) and disorder thus modifies the physics at a microscopic level. We analyze 
the coherence properties of the cloud through measurements of the momentum distributions, for 
two disorder strengths, as a function of its degeneracy. For moderate disorder, the emergence of 
coherence remains steep but is shifted to a lower entropy. In contrast, for strong disorder, the growth 
of coherence is hindered. Our study is an experimental realization of the dirty boson problem in a 
well controlled atomic system suitable for quantitative analysis. 
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Together with band structure and interactions, disor- 
der is a key ingredient for the understanding of transport 
in condensed matter physics At low temperature, it 
affects the conductivity of a metal and it even induces 
phase transitions to insulating states 0- A striking ex- 
ample is Anderson localization 0], which has recently 
been observed in 3D ultracold gases [3l- 

Disorder is especially relevant in 2D systems, such as 
Si-MOSFET [il, or thin metal films 0, in which quan- 
tum phase transitions to insulating phases have been 
observed. Moreover, in high-Tc superconductors, dop- 
ing intrinsically introduces inhomogeneitics in the CuO- 
planes • Understanding the complex interplay between 
disorder and interactions in these systems remains a ma- 
jor challenge. 

Whereas the above mentioned electronic systems are 
fermionic, superconductivity originates from the bosonic 
nature of Cooper pairs. As long as disorder docs not 
break the Cooper pairs, the problem is reduced to a 
study of dirty bosons [1, Q . It has mainly been stud- 
ied numerically in the framework of the disordered 2D 
Bose-Hubbard model. Disorder can both favor or disfa- 
vor superfluidity [l^l, and the occurrence of a Bose glass, 
an insulating, gapless, compressible phase has been pre- 
dicted Hi. 

In the context of ultra-cold atoms, the properties of 
disordered trapped Bose gases have been studied both 
in ID and 3D [Ij. In 2D Bose gases in the 

absence of disorder, the Berezinskii-Kosterlitz-Thouless 
(BKT) superfluid phase transition [l^ has been experi- 
mentally studied through the modification of the gas co- 
herence prop erties associated with the pairing of thermal 
vortices [3 ■ In a continuous system, the effect of 
disorder on the BKT superfluid transition is expected to 
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depend on the correlation length of the disorder a, which 
has to be compared to the characteristic length scales of 
the cloud such as the thermal de Broglie wavelerigth AdB 
and the healing length ^, i.e. the vortex core size [l8|. For 
slowly varying disorder (cr ^ ^, Ads), the physics can be 
locally described by the homogeneous BKT transition, 
and disorder causes a percolation transition of superfluid 
islands. In contrast, for a microscopically correlated dis- 
order (ct < ^, Adfi), tunneling is possible and the very 
nature of the phase transition is affected. 

A recent Monte- Carlo study of the homogeneous 3D 
Bose gas in the presence of speckle disorder [l^ has 
shown that, depending on its correlation length a, the 
disorder can cither reduce the critical temperature (due 
to quantum localization) or increase it because of the 
reduction of the available volume (see Fig. 14 in [l^ ) . 
However there have been so far no theoretical prediction 
for the effect of disorder on the BKT superfluid tran- 
sition in a continuous 2D system. As in 3D, Anderson 
localization Q and percolation phenomena are likely to 
affect the superfluid transition, but reducing dimension- 
ality should enhance their effects [l^, [l^ . In addition, 
new phenomena affecting specifically the BKT transition 
such as enhanced phase fluctuations or vortex pinning 
j20| may play an important role. 

In this letter, we present an experimental study of 
the effect of microscopically correlated disorder on the 
coherence properties of a 2D ultracold atomic gas near 
the BKT superfluid transition. As in [2l|, the coherence 
properties are probed by the study of the momentum dis- 
tribution, which is the Fourier transform of the first order 
correlation function gi [22| ■ We observe that an adiabatic 
ramping up of the disorder results in a suppression of the 
low momentum peak, i.e. a decrease of coherence. In par- 
ticular, for a moderate disorder strength of 0.4 times the 
temperature, we measure a small shift of the emergence 
of coherence towards low entropy. For stronger disorder 
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strength of the order of the temperature, the growth of 
coherence is significantly hindered both as a function of 
entropy and temperature. 

Our experiment starts with a quantum degenerate 2D 
Bose gas in a trap obtained from a combination of a blue 
detuned TEMqi beam, which confines the gas in a hori- 
zontal plane, and a red detuned Gaussian beam for the 
in-plane confinement [2l|. The trap oscillation frequen- 
cies are lOx/'^t^ = 8 Hz, ujy /2ti = 15 Hz, lOz/Itt = 1.5 kHz. 
The atom number N is varied between 2 x lO"' and 
6 X 10'^ in order to change the degeneracy of the gas 
across the BKT transition. The temperature, measured 
from a fit to the wings of the momentum distribution us- 
ing a Hartree-Fock mean- field model [2l|, remains con- 
stant at 64.5±2.0nK. At this temperature, ~ 70% of the 
atoms are in the ground state of the vertical harmonic 
oscillator. The dimensionless 2D interaction strength is 
g = y/STras/a^ = 0.096, where Og = 5.3 nm is the 3D 
scattering length, = \/ hjmujz ~ 0.28 /im the verti- 
cal harmonic oscillator characteristic length, m the atom 
mass, and h the reduced Planck constant. 

The disorder potential is a speckle pattern produced 
by a 532 nm laser beam, which passes through a diffusive 
plate and is focused on the atoms. The repulsive disor- 
der potential is characterized by its mean value V (equal 
to its standard deviation) and its correlation lengths, 
inversely proportional to the numerical aperture of the 
optical system [23|- Given the intensity of the beam 
and its transverse waist radius of 1 mm, the maximum 
value of V felt by the atoms is Fmax = fee x 60(10) nK. 
As the beam is tilted by 30°, the in-plane disorder is 
effectively anisotropic [2J]. The correlation lengths of 
the disorder are such that (Tx/2 = cry = 0.5 /^m (half- 
width at 1/ y/e) . These correlation lengths are of the or- 
der of both the thermal de Broglie wavelength AdB = 
\j2'iTh'^ /mk^T sa 0.73 yum and the healing length (at 
the BKT transition) ^ = Xdn/ V^cg ~ 0.82 /im, where 
Dc « log(380.3/g) w 8.3 is the BKT critical phase space 
density |25| . However, (Tx and ay are small compared to 
the Thomas- Fcmn radii of the cloud at the BKT tran- 
sition Ix = |^^^=25/xm, ly = {uJx/ujy)lx = 13 /im, a 
necessary condition for self-averaging measurements. In 
our experiment, we use a single realization of the speckle 
pattern. 

In the experimental sequence, the disorder potential 
is slowly ramped up in 250 ms after the preparation of 
the 2D gas. After a holding time of 250 ms, all trapping 
potentials including the disorder are switched off and the 
atom cloud expands in 3D during a free fall of 83.5 ms. 
The column density of the gas along z is then measured 
by fluorescence imaging from the top. As explained in 



21| , it reflects the in-trap momentum distribution in the 
x, y plane. Since, the momentum distributions appear 
to be cylindrically symmetric, we perform an azimuthal 
averaging [26j to obtain the momentum profiles n{k) as 
a function of the wavenumberfc. 

We study the effect of disorder on the momentum dis- 
tribution for different initial conditions, both above and 
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FIG. 1. (Color online) Azimuthally averaged momentum dis- 
tribution profiles for iV = 5.6 x 10* (a), iV = 3.8 x lO'* (b) 
and N = 2.9x lO" (c). In each case, we show the influence of 
the disorder: no disorder V = (black triangles), OAVmmx 
(green open squares), V^max (red circles). 



below the BKT transition. A narrow momentum distri- 
bution is the signature of a large coherence length in 
the sample. We quantify the degeneracy of the non- 
disordered gas with the ratio N/Nc, where Nc is the 
critical atom number of an ideal 3D Bose gas in our 
anisotropic trap. From our previous study of the 2D 
Bose gas in the absence of disorder [2l|, we know that 
the BKT phase transition happens at N/Nc « 1.26 for 
our parameters. In Fig.[TJ we compare the momentum 
distribution without disorder to the results obtained af- 
ter ramping up the disorder potential to V — 0.4 I^max 
and V Fmax- For N/N^ 1.06 (Fig.[lJ;) and in the ab- 
sence of disorder, the gas is in the normal phase. In this 
case, the addition of the disorder has little effect, reduc- 
ing slightly the low momentum population (fc < 2 /im~^). 
For N/Nc = 1.32 (Fig.[T})), the gas has just entered the 
superfluid phase in the absence of disorder and a low mo- 
mentum peak is clearly present. In this case, the disorder 
has a strong effect. The population at very low momen- 
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turn (fc < 0.5 fj,m.~^) is strongly suppressed and the low 
momentum peak almost disappears. The profiles with 
disorder are then qualitatively similar to the one in the 
normal phase (Fig.[ll:). For N/N^ = 1.75 (Fig.[I]a) and in 
the absence of disorder, the gas is deep in the superfluid 
phase with a large low momentum peak. In this case, 
the addition of the disorder leads to a reduction of the 
height of the peak but not to its disappearance. In all 
our data, adding disorder always results in a reduction 
of the coherence of the Bose gas. 

It should be noted that applying slowly the disorder 
preserves the entropy. When the disorder potential is 
ramped up to a mean value l^max in 250 ms and then 
down in 250 ms, we find no heating and no atom loss com- 
pared to the non-disordered situation, within our exper- 
imental precision (±lnK). This uncertainty comes from 
both the shot-to-shot fluctuation of the experiment and 
the accuracy of the fitting procedure. We thus observe 
that adding the disorder is a reversible process. In the 
following, we assume that consecutive pictures with and 
without disorder correspond to the same entropy (Fig.[T|). 
Presenting our data as a function of entropy is then a 
natural choice. 

Even in the absence of disorder, we however do not 
have experimental access to the entropy. To find the cor- 
respondence between the ratio TV/TVc, extracted from our 
measurements and entropy, we rely on quantum Monte- 
Carlo simulations of the non-disordered in-situ distri- 
bution (27j . from which the entropy can be determined 
because of the scale invariance of the 2D Bose gas [28j . 
The calibration of the average entropy per particle S as 
a function of N/Nc for our experimental conditions is 
shown in the inset of Fig.O 

In order to analyze our result in a simple way, we would 
like to characterize the degree of coherence of the gas with 
a single number and not with the full momentum distri- 
bution. A natural quantity to consider is the coherence 
length, i.e. the inverse of the width of the momentum 
distribution. However, in our case, the width of the mo- 
mentum distribution saturates because of the limited res- 
olution of our imaging system for highly coherent clouds 
[2l| . As an alternative, we choose to focus our analysis 
on the fraction of atoms Nq/N in the central pixel of the 
momentum distribution (fc < 0.2 /im~^). It is also a well- 
defined model-independent quantity and it is related to 
the fraction of atoms that are coherent on a length scale 
larger than ^5^m [2^. We plot Nq/N as a function of 
the entropy per particle S (Fig. [2]) and find that, at fixed 
entropy, the coherence of the gas is reduced in the pres- 
ence of disorder. For V = 0.41/max, the emergence of 
coherence is slightly shifted to a lower entropy per par- 
ticle, by 0.2(1) ks, compared to the non-disordered case. 
For V = Vmaxj we never reach a sufficiently low entropy 
to observe a large increase of coherence. 

Since the phase diagram of disordered systems is typi- 
cally presented as a function of disorder and temperature 
[soj . we now complement our analysis of the coherence 
as a function of these quantities. This means that we 
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FIG. 2. (Color online) Fraction of atoms Nq/N in the central 
pixel of the momentum distribution as a function of the aver- 
age entropy per particle S: non-disordered case V = (black 
triangles), V = 0.4Fniax (green open squares), V — Vmax 
(red circles). Each point results from the averaging of 5 ex- 
perimental profiles and the error bars are statistical. The line 
corres pon ds to a Monte-Carlo simulation in the absence of dis- 
order [22] . Inset : entropy per particle measured by quantum 
Monte-Carlo simulations at 64.5nK (blue cross) and fitted by 
a 3^'^ order polynomial (black line). The dashed lines indicate 
the BKT transition. 



have to determine the temperature from the experimen- 
tal disordered profiles. In the absence of an exact the- 
oretical model for 2D disordered gases, we use our non- 
disordered Hartrce-Fock mean-field model ^l| , which we 
expect to be valid at large momenta. Experimentally, 
we fit to the wings of the distribution between a variable 
cut-off momentum kc and 12/im~^. We find that, for 
2.75 /im^^ < kc < 3.75 ^m~^, our signal to noise ratio 
is sufficient and the fitted temperature varies typically 
less than 1 nK, indicating that our model is reasonably 
accurate in this range. We use kc = 3.5/im~^ in the 
following analysis. For V = l^max the temperature is 
found to increase on average by 5.5 nK compared to the 
non-disordered case (T ~ 64.5 nK). 

Figure[3] presents Nq /N as a function of the temper- 
ature normalized to Tc, the critical temperature for an 
ideal 3D Bose gas in our anisotropic trap and the mea- 
sured atom number. The results without disorder and 
with a disorder of amplitude 0.4 l^max sltc similar and no 
clear shift is visible. Within our accuracy, we can con- 
clude that for this amount of disorder the coherence prop- 
erties are weakly affected. Note that the classical perco- 
lation threshold ^ 0.52 V ~ 13 nK [3l| is larger than the 
mean-field chemical potential at the phase transition in 
the absence of disorder ^c — gf^^Dc/mXi^'^ w fee x 8nK. 
The observation of a weak effect of disorder is a signature 
that classical percolation is not relevant in the regime of 
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FIG. 3. (Color online) Fraction of atoms Nq/N in the central 
pixel of the momentum distribution as a function of the nor- 
malized estimated temperature for disorder strengths, V — 
(black triangles), V = OAVmax (green open squares), and 
V — Vmax(red circles). is the critical temperature of an 
ideal 3D Bose gas in our anisotropic trap. Each point results 
from the averaging of 5 experimental profiles and the error 
bars are statistical. The line corresponds to a Monte-Carlo 
simulation in the absence of disorder [2j . 



a microscopically correlated disorder (cr ^, Ade)- 

For V = T^max, the coherence increases much slower 
when T/Tc decreases. The coherence properties of the 
gas are greatly modified. This finding contrasts with the 
3D quantum Monte-Carlo calculation, which predicts a 
significant effect of the disorder only at larger disorder 
strength We thus show the enhanced role of dis- 

order in 2D as compared to 3D. At our lowest temper- 
atures, Nq/N does not reach the value of 0.0035 which 
corresponds to the superfluid transition in the absence 



of disorder [2T|. Although we do not directly measure 
the superfluid fraction [ITj . since the superfluid transi- 
tion is generally associated with the apparition of long 
range coherence, we can suspect that our system is not 
in a superfluid phase and that the critical temperature 
for superfluidity is shifted down by a sig nificant amount. 
Actually, for our amount of disorder [32| , the existence of 
a superfluid is not guaranteed even at zero temperature 
because of the disorder-driven quantum phase transition 
from a superfluid to an insulating Bose glass phase [soj . 

In conclusion, we have shown that a microscopically 
correlated disorder (cr ^ Ci-^dB) always reduces the co- 
herence of a 2D Bose gas, both at constant temperature 
and entropy. For moderate disorder strength of 0.4 T, 
the reduction is weak although we are able to measure 
a small shift of the emergence of coherence toward low 
entropy. A disorder strength of the order of gas temper- 
ature leads to a qualitative change of behavior with a 
suppressed coherence growth. We interpret the observed 
strong suppression of the coherence growth as a large 
shift of the superfluid transition. Theoretical studies in 
our experimental conditions would be of great interest 
and would permit to strengthen our analysis. The mech- 
anism of the disorder action can also be addressed both 
experimentally and theoretically. In the future, similar 
studies on strongly interacting fermions would add the 
possibility to elucidate the physics of disorder-induced 
breaking of bosonic pair. 
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N.B.: In the last stage of the redaction of our paper, 
we have learned about a complementary work about dis- 
ordered 2D Bose gases in the deep superfluid regime [s^ . 
It is focused on the different behaviors between quasi- 
condensate fraction and coherence. 
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